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ABSTRACT 

It is long debated if pre-biotic molecules are indeed present in the interstellar medium. De¬ 
spite substantial works pointing to their existence, pre-biotic molecules are yet to be discovered 
with a complete confidence. In this paper, our main aim is to study the chemical evolution of 
interstellar adenine under various circumstances. We prepare a large gas-grain chemic al network 
by considering various pathways for the formation of adenine. Maiumdar et al. ( 20I2ll proposed 
that in the absence of adenine d etection, on e could try to trace two precursors of adenine, namely, 
HCCN and NH 2 CN. Recently iMerz et al.l ( 20I4ll . proposed another route for the formation of 
adenine in interstellar condition. They proposed two more precursor molecules. But it was not 
verified by any accurate gas-grain chemical model. Neither was it known if the production rate 
would be high or low. Our paper fills this important gap. We include this new pathways to find 
that the contribution through this pathways for the formation of Adenine is the most dominant 
one in the context of interstellar medium. We propose that observers may look for the two pre¬ 
cursors (C 3 NH and HNCNH) in the interstellar media which are equally important for predicting 
abundances of adenine. We perform quantum chemical calculations to find out spectral proper¬ 
ties of adenine and its two new precursor molecules in infrared, ultraviolet and sub-millimeter 
region. Our present study would be useful for predicting abundance of adenine. 


Subject headings: Bio-molecules, Interstellar Medium, Astrochemistry, Molecular cloud, Star formation 


1. Introduction 


According to the Cologne Database fo r Molec¬ 


ular Spectroscopy ( Muller et al. 2001 . 2005[l . 


around 180 molecules have been detected in the in¬ 


terstellar medium (ISM) or in circumstellar shells. 
Among these, severals species are organic in na¬ 
ture. Existence of these complex molecules could 
not be explained wit hout a proper consideration 
of interstellar dusts (IHasegawa, Herbst fc Leung 
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1992; Das . Acharvva fc Chakrabartil20f^ Booeert fc EhretrfiteiiDdl rform a comparative study among various 


2004 : lGibb et al.l 2004li . Several attempts 


were 


made over the past few years to study physical 
and chemical processes on the interstellar grains 


2007 

Das et al. 

201(1 

Das et al. 


2008alpas, Acharvva fc Chakrabarti 
2013bt1sahuff2015 b Several exper¬ 


imental results have been reported in explaining 
the importance of various surface pr o cesses on in¬ 
terste llar dusts ( loPDolo et al.l 2008t lOberg et al.l 
2 OO 9 I) . making it imperative to incorporate the 


grain surface chemistry extensively while predict¬ 
ing their abundances. 

The problem of Origin of life is a long standing 
puzzle and formation of ami no acids in th e labora¬ 
tory by well known work of Milleil ( 1953ll ushered 
a new direction of research in this area. Mo re re¬ 
cently, IChakrabarti fc Chakrabaidil (j2000al lt3l for 
the first time, suggested that perhaps the process 
of formation of the complex molecules such as 
Adenine and other constituents of DNA is very 
generic and such complex pre-biotic molecules 
should be formed during any star forming pro¬ 
cess. I n the absence of proper reac t ion cro ss sec¬ 
tions, IChakrabarti fc Chakrabai^ ( 2000 atl used 
neutral-neutral reaction rates to compute ade¬ 
nine abundance with successive addition of HCN. 
This was later impr oved upon with more real¬ 
istic cross sections (jChakrabarti fc Chakrabarti 


2000 bh and more realistic abundance was ob¬ 
tained. Presence of interstellar grains were in 
explicitly added in these works, and were incorpo¬ 
rated indirectly by a higher rate coefficients in 


forma tion. A follow up study bv iMaiumdar et al 


( 20121 1 explicitly considered presence of grains and 
showed that this significantly alters the abun¬ 
dance of adenine. They carried out several pre¬ 
scriptions of rate coefficients for successive HCN 
addition reactions. They also considered radical- 
radical/radical-mol ecular reactions proposed by 
Cupta et al. ( 201lh for the formation of adenine. 


Their results suggest that radical-radical/ radical- 
molecular reactions dominate s over the neutral - 
neutral reactions. Recently iMerz etldl ( 2014 1 
used the concept of retro-synthetic analysis to pro¬ 
duce interstellar adenine from observed interstel¬ 
lar molecules such as C 3 NH, HNCNH and its iso¬ 
mer H 2 NCN. They used MP2/6-311-I—|-C(2d,2p) 
method to calculate various chemical parameters 
involving a six step mechanism. This motivated us 


available pathways for the formation of adenine in 
interstellar region. This we present below. 

The plan of this paper is the following. In Sec¬ 
tion 2, computational methods are discussed. Re¬ 
sults are presented in Section 3, and finally in Sec¬ 
tion 4, we draw our conclusions. 


2 . 


2 . 1 . 


Methods and Computational Details 
Chemical modeling 


We develop a large gas-grain chemical network 
to explore chemical evolution of interstellar ade¬ 
nine and its related species. For the gas-phase 
chemical network, we follow UMIST 2006 data 
base. Formation of adenine via various reac¬ 
tion pa thw ays used inlc hakrab arti fc Chakrabarti 
( 200(1 al fhll , Cunta et al.l (|2011 ) and Merz et al. 
(120141 ) are included. Our gas phase chemical net¬ 
work consists of more than 6000 reactions between 
650 species. 


low 

Haseffawa. Herbst & Leuna 

(1992): 

Cuopen & HerbstI 

(20oi; 

(HI; 

Das. Acharvva & Chakrabartil ( 

20101): Das & Chakrabarti 

Das et al. (2013b. 

20151) 

. Our surface reac- 


tion network contains 292 reactions. We consider 
gas-grain interaction to frame the actual inter¬ 
stellar scenario. Cas phase species are allowed to 
accrete on interstellar grains. Crains are actually 
acting as a catalyst for the formation of several 
complex interstellar species. Formation of the sim¬ 
ple and most abundant interstellar molecule, H 2 
cannot not be explained without the consideration 
of interstellar grain chemistry. Binding energies 
of surface species mainly dictate chemical com¬ 
position of interstellar grain mantle. We consider 
the most upd ated intera c tion ba rrier energies as 
mentioned in IPas et a~ ( 2013lJ) and references 
therein. Depending on the energy barriers, sur¬ 
face species would move throughout the grain sur¬ 
face by thermal hopping or tunneling whichever 
is faster. At low temperatures, for lighter species 
such as H atom, tunneling is much faster. During 
these movements, surface species would react with 
any suitable reactant partners. Based on their en- 


evaporate ( 

Haseeawa. Herbst & LeunaIl992l) and 

populate the 

gas phase. Cosmic rav induced 

evaporation ( 

haseaawa & Herbst 

1993) and non- 

thermal desorption ( 

Carrod. Wakelam & Herbst 
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Table 1: Available and estimated rate co-efficients for the formation of interstellar adenine in the gas phase 


via various reaction pathways 


Pathways used in 

Reactions 

Rate coefficients 


(i)HCN + HCN ^ CH(NH)CN 

8.38 X 10“^*^cm^ sec“^ 

Chakrabarti Chakrabarti f2000a.bb 

(m)CH(NH)CN + HCN ^ NH 2 CH(CN )2 

3.43 X 10“^^ cm^ sec“^ 

Majumdar et al. f2012^ 

(iii)NH 2 CH(CN )2 -1- HCN -)■ NH 2 (CN)C = C(CN)NH 2 

3.30 X cm^ sec""^ 


(iv)NH 2 (CN)C = C(CN)NH 2 -I- HCN -)• C 5 H 5 N 5 

3.99 X sec""^ 


(v)HCCN -(- HCN -> Molecule 1 

2.13 X CHI'* sec“^ 


(vi)Molecule 1 + H —>• Molecule 2 

7.96 X 10~^ cm^ sec“^ 


(vii)Molecule 2 + NH 2 CN —Molecule 3 

6.24 X 10”^^ cm^ sec“^ 

^^^t^^t al. f 2011 ) 

(viii) Molecule 3 + CN Molecule 4 

1.80 X 10“®cm^ sec“^ 


(ix)Molecule 4 + H ^ Molecule 5 

8.71 X 10“®cm^ sec~^ 


(x)Molecule 5 -I- CN -> C 5 H 5 N 5 + HNC 

1.89 X 10“®cm^ sec“^ 


(xi)Molecule 5 -I- CN ->■ C 5 H 5 N 5 -I- HCN 

1.91 X 10“^cm^ sec“^ 


(xiqCsNH + HNCNH ^ C 4 N 3 H 3 

8.30 X 10“^*^cm'* sec“^ 


(xiii)C 4 N 3 H 3 -1- HNCNH ->■ C 4 H 3 N 3 -I- NH 2 CN 

6.43 X sec“'^ 

Merz_et_^ 1 ^^ 01 ^ 

(xiv)C 4 H 3 N 3 + HNCNH C 5 N 5 H 5 

1.36 X 10“*’sec“i 


(xv)C 5 N 5 H 5 -I- HNCNH -s- N 5 C 5 H 5 + NH 2 CN 

1.00 X 10“®cm^ sec“^ 


(xvpNsCsHs -1- HNCNH C 5 H 5 N 5 -1- NH 2 CN 

1.89 X 10“^^sec“^ 
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Fig. 1.- Comp arison of results of v arious adeni n e form a- 

tion pathways. Pathways proposed by IMerz et al.l I I2014D is 
dominating over all other pathways. 


2007 : Das et al. 20151) mechanisms are also con¬ 
sidered in our model. So, in brief, gas and grains 
are interacting with each other to exchange their 
chemical components by various means. Since all 
the processes are random, Monte Carlo method 
would be appropriate to use while dealing with 
this randomness. However, it requires huge com¬ 
putational time to compute with our vast gas 
phase and grain phase chemical networks simulta¬ 
neously. Thus, we use traditional Rate equation 
method to handle our chemical network on grains. 
Detailed discussions on our gas-grain chemical 


mode l are already presented in [D as et al.l (l2013b 


2015 ): Maiumdar. Das fc Chakrabartil ( 2014allb ). 


All the reactions, which are considered here for 
the formation of adenine are shown in Table 1. 
Maiumdar et al. ( 2012 . 201, Sh considered neutral- 
neutral pa thways (reaction number (i)- (iv) of Ta- 


of 


ble 1) 

and radical-radical / radical- molecular 


Chakrabarti fc Chakrabartil ( 2000al lb) 

pa t hways 


(reaction numbers (v)-(xi)) of lCunta et al. ( 2011 ) 
and concluded that the production of adenine 
is dominated by radical-radical/radical-molecular 
reaction pathways. They identified HCCN and 
NH 2 CN as two main precursor molecules, which 
are responsible for the production of adenine. 
HCCN is highly abundant in the interstellar 


space (j.Tiurvj 
Formation of 


■ig nly a 

J l200(it 


Guelin fc Cernicharol [19911). 


HCCN on the interstellar grains 


were al so considered ((Hasegawa . H erbst fc Leung 
I 992 I) . McGonagale fc Iryind ( 1996 ) conducted a 


deep search for HCCN towards TMC-1 and sey- 
eral CMC’s yia its N(J) = 1(2) 0(1) transition. 
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They came up with an upper limit of fractional 
abundance with respect to molecular hydrogen to 
be ^ 2 X 10“^°. Existence of NH 2 CN in the in¬ 
terste llar cloud was first reported bv I Turner et al.l 
( I 975 II . Subsequently, it w as observed in both di f¬ 
fuse and dense cloud s by Liszt fc Lucas I ( 2001 1. 
Woodall et al.l ( 2007l l predicted a steady state 
fractional abundance of 2.02 x 10“^° for NH 2 CN 
with respect to H 2 . 

Merz et ^ (12014 1 proposed a new pathway (re¬ 


action nos. x(ii)-(xvi) of Table 1) for the pro¬ 
duction of adenine. They used retro synthetic 
analysis by using two new species (C 3 NH and 
HNCNH). Carbodimide ( HNCNH), is aii i somer 
of cyanamide (NH 2 CN). iMcGuire er al.l ( 2013 1 
obtained the abundance of this molecule from ice 
mantle experiments. They proposed that tau- 
tomerization of NH 2 CN on dust grain ice mantles 
is the dominant formation pathway for HNCNH. 
They found its abundance would be 10% of 
NH 2 CN (in Sgr B2(N) column density of NH 2 CN 
is ^ 2 X 10^^ cm“^). Since, this was below the 
detection limit of any current astronomical facil¬ 
ity, they proposed to observe HNCNH by those 
transitions which are amplified by masing. In our 
chemical model, we consider that 10% of NH 2 CN 
could be converted into HNCNH. 

We us e semi -empirical relationship developed 
bv iBatesI ( 1983ll for the computation of rate co¬ 
efficients of the exothermic and barrier less reac- 


tion p athways (reaction nos. xii-xiv) o n Merz et al 

(I 2 OI 4 II . 


K = lxlO"^^H^(6Eo-hA^-2)(3"-^V(3A^-7)!cTO^s"^ 

( 1 ) 

where, Eq is the association energy in eH, is 
the transition probability (in s“^) of the stabiliz¬ 
ing transition (the numerical value of which may 
be taken to be 100 unless better information is 
available) and N is the number of nuclei in reac¬ 
tants. For the reaction numbers xii, xiii & xiv, 
we use association energies of —5.2, —0.6, —72.8 
kcal/mol respectively. 

If the calculated rate coefficients from Eqn. 1 
exceeds the limit set by the following equation 
(Eqn. 2), then this limiting value was adopted: 


K = 7.41 X 10”^V/^ 


( 10 //.) 




( 2 ) 


where, a is the polarizability in A°^, is the re¬ 
duced mass of the reactants in amu scale as 


suggested bv iBated ( 1983 1. 

Rate coefficients of reactions having activation 
barriers (reaction nos. xv-xvi) are calculated by 
using conventional transition state theory. Ac¬ 
cording to this theory, the rate coefficient has the 
following form: 


k{T) = (KBT/hCo)exp{-AG/RT) s“ 


(3) 


where, Kb is the Boltzmann constant, h is the 
Plank’s constant, T is the temperat ure, Cn is the 
concen tration (set to 1 , by following I Jalbout et al 


(|2008f) l. R is the ideal gas constant, and AG is 
the free energy of activ ation. From the q uantum 
chemical calculations bv iMerz et al. ( 2014ll . we are 
having AG = -1-0.4 kcal/mol and -1-1.1 kcal/mol 
respectively for reaction numbers xv and xvi. 

2.2. Spectroscopic Modeling 

An educated estimate of spectral properties is 
essential before observing any unidentified species. 
It is reported in earlier astrophysical literature 
that DFT and TDDFT (jPunge fc Cross! Il984ll 
computation methods would be applied to var¬ 


ious astrophysica l p roblems (Puletti et al. 2010 t 


Pieve et al. 2014 1. Das et al] ~( 2015h discussed 


the necessity of quantum chemical calculations 
prior to any spectroscopical survey. Since our 
chemical model would predict the abundances 
of adenine and its related species, we believe it 
would be useful to present various spectral as¬ 
pects of these species. For finding the various 
spectral aspects, we perform quantum chemical 
calculati ons by using Caussian 09W program. B y 
following lMaiumdar. Das fc Chakrabartil ( 2014bl l , 
for the vibrational frequencies, we use B3LYP/6- 
311-1—|-C(d,p) level of theory. In case of grain 
phase species, Polarizable Continuum Model 
(PCM) is used with the integral equation for¬ 
malism variant (lEFPCM) as the default SCRF 
method. For the electronic absorption spectrum, 
we use time dependent density functional theory 
(TDDFT study). Since most of the gas phase com¬ 
plex molecules were observed by their rotational 
transitions in the mm or sub-mm regime, we carry 
out quantum chemical calculations to find out the 
rotational transitions of our desired species. Com¬ 
putation of the anharmonic frequencies require the 
use of analytic second derivative of energies at the 
displaced geometries. But the CCSD method in 
Caussian 09W program only implements energies, 
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Table 2: Initial elemental abundances 


Species 

Abundance 

H 2 

5.00 X 

He 

1.00 X 10“°'^ 

N 

2.14 X lO"®*^ 

0 

1.76 X 10“°^ 

H 3 + 

1.00 X 10"“ 

c+ 

7.30 X 10"°'^ 

s+ 

8.00 X 10"°® 

Si+ 

8.00 X 10"°° 

Fe+ 

3.00 X 10"°° 

Na+ 

2.00 X 10"°° 

Mg+ 

7.00 X 10"°° 

p + 

3.00 X 10"°° 

C1+ 

4.00 X 10"°° 

e~ 

7.31 X 10"°° 

HD 

1.6 X 10"°° 


so analytic second derivative of energies are not 
available at this level of theory. So there are no 
options in Gaussian 09W program to compute the 
rotational and distortional constants at the CCSD 
level of theory. Here, we use B3LYP/aug-cc-pvTZ 
level of theory, which are a l so proven to be v ery ac¬ 
curate ( Carles et al. 2013 : IPas et al.llioisll . This 
level of theory would be appropriate because there 


are so me earlier studies (such as, iBrunken et al 


( 2006ll ) for finding rotational and distortional con¬ 
stants of Uracil (one nucleo base of RNA beside 
adenine, cytosine and guanine). Corrections for 
the interaction between the rotational motion and 
vibrational motion along with the corrections for 
vibrational averaging and an-harmonic corrections 
to the vibrational motion are also considered in 
our calculations. Obtained rotational and distor- 


tional constants are then used in SPCAT (Pickett 


19911 1 program to predict various rotational tran¬ 
sitions. Output of the SPCA T program is then 


direc tly used in ASCP Prgram ( Kisiel et al.lll998L 
I 2 OOOII to find out the rotational stick diagram of 
the desired species. 


3. Results and Discussion 

3.1. Results of Chemical modeling 



Chakrabarti fc Chakrabaiui ( 2000al f9l : Maiumda 


(|2ni2[ l considered successive HCN addit i on for 


the formation of adenine. Cupta et al.l ( 2011 1 


proposed adenine form ation pathways sta rting 
with HCCN and HCN. Gupta et al.l ( 201 ill also 
pointed out that their reaction pathways would 
also be feasible in the grain phase. A com¬ 
pletely new pathway has been proposed recently 


Fig. 2.- Radial distribution of density profile. 
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by iMerz et alJ (l2014l) . where adenine could be 
produced without any HCN addition. So, we 
wish to compare four different pathw ays where the 


first o ne is the gas phase pathway of lGuota et al 


( 201 lh . t he second one is the pathw a y cons id- 


ered bv Ghakrabarti & Ghakrabarti 

(l 2000 allb) 

and updated bv Maiumdar et al. ( 

2012 h. the 


thi rd one is the ga s phase pathway proposed 
by Merz et al. ( 20141 ). As a fourth method, we 
consider both gas and grain phase pathways of 
iGupta et al. ( 2011 ). Since, no mention wa s made 
of t he grain phase reactions in Merz et al.l ( 2014h 
and Chakrabarti &: Chakrabartil (|2000ai lJ). we do 
not consider these pathways in the grain phase. 
In general, we consider that the gas and the 
grain interact with each other to exchange their 
chemical components. So, grain phase could be 
populated by the accretion from the gas phase 
and gas phase could be populated by the grain 
phase species through evaporation mechanisms 
such as, thermal desorption, non-thermal des¬ 
orption and cosmic-ray induced desorption. Ini¬ 
tial elemental abundances with respect to to¬ 
tal hydrogen nuclei (Tab le 2) are consider ed to 
be similar to what is in Das et al.l (gOldb ) and 


Maiumdar. Das fc Chakrabartil ( 2014bl) . These 
are the typical low metallic abundances which are 
often adopted for TMC-1 cloud. In Fig. 1, a com¬ 
parative study between these four pathways are 
shown for the chemical evolution of adenine. To 
draw this, we take number density of hydrogen nu¬ 
clei existing in all possible forms {uh) to be = 10 "'^ 
cm“^. Ay = 10 and T = 10 K. Various evap¬ 
oration mechanisms are considered here, among 
them, around low te mperatures, non-thermal des¬ 
orption mechanism (IGarrod, Wakelam fc Herbst 


2007t iDas etZll2015[) is the most efficient means 
to transfer these complex molecules to the gas 
phase. A moderate value of the desorption 
parameter 0.05) is considered. Significant 
difference in the abundance of adenine could 
be observed betwe en the two considerations of 
iGupta et al. ( 2011). In on e case, only gas phase 


pathways of Gupta et al. ( 201ll) is considered 


and in another case, the same pathways are 
considered for both the phases (gas and grain). 
From Fig. 1, we find that the peak abun¬ 
dances of adenine for these two cases appear 
to be 5.3 X 10“^^ and 1.5 x 10“^^ respectively. 
Since grain phase production is efficient, produc¬ 


tion is significantly enhanced in the case when 
grain phase pathway s are also considered. While 
using pathways of iGhakrabarti fc Ghakrabai^ 
( 2000 al fa). peak abundance of adenine turns 
out to be 4.8 x 10“^^ and u sing pathways pro- 
(|2014l ). the peak abun- 
10“^"^. Using neutral- 


posed by iMerz et al. 
dance becomes 2.5 x 


neutral rate coefficients in a n isoth ermal cloud, 
Ghakrabarti fc Ghakrabarti ( 2000al ) found the 
abundance of ~ 5 x 10“^^, one magnitude higher 
than this. It is also clear fro m Fig. 1 that the 


adenine formation pathways of I Merz et al.l ( 2014h 
dominates over all other pathways. In Fig. 1, we 
show chemical evolutions of HNCNH and C 3 NH 
(which are required for the f ormation of adeni ne 
in the pathways proposed bv iMerz et al.l (|2014l) '). 
These molecules could be used as the precursors 
for estimating the abundances of adenine in inter¬ 
stellar region. 

3.1.1. Mixing model 

Chemical complexity of any interstellar region 
is highly dependent on surrounding physical con¬ 
dition. Formation of bio-molecules would also be 
influenced by surrounding physical process. To 
describe a realistic situation, we prepare a spe¬ 
cial model named “Mixing model”. We consider 
that density pr ofile f ollows p ^ distribution 


as described by Shul ( 1977t) . At the outer bound¬ 
ary of the cloud the density is 10 ^ cm“^, inner 
boundary location is chosen in such way that the 
density would become 10® cm“® (Fig. 2). For 
concreteness, we choose the outer boundary to be 
at 15962 AU and the inner boundary becomes at 
159.62 AU. To reduce computational time we fur¬ 
ther subdivide this cloud into four shells, inner¬ 
most shell is extended from 159.62 AU to 504.76 
AU which has an average density of 2.12 x 10^ 
cm“®. The shell no. 2 has an average density of 
2.12 X 10® cm“® and extends from 504.76 AU to 
1596.2 AU. The shell no. 3 has an average den¬ 
sity of 2.12 X 10® cm“® and extends from 1596.2 
AU to 5047.63 AU. The fourth shell extends from 
5047.63 AU to 15962 AU and has an average den¬ 
sity of 2.12 X 10'’’ cm“®. Since the cloud is dense, 
we assume visual extinction (Ay) of 10 and tem¬ 
perature 10 K throughout the cloud. 

We assume that a systematic mixing is going 
on throughout the cloud. Matter of shell 4 would 
contribute {p 4 % of shell 4 matter) to shell 3, ma- 
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terial of shell 3 would contribute (p 3 % of shell 3 
matter) to shell 2 , material of shell 2 would con¬ 
tribute {p 2 % of shell 2 matter) to shell 1 and due 
to re-entry of the outflow at the outer edge, some 
matter would re-enter {pi% of shell 1 matter) into 
shell 4. Each transport should have different time 
scales. Inward material transfer (from shell 4^ 
shell 3—shell 2 —>-shell 1) is assumed to take place 
with the sound speed of /nip ^ 2.873 x 10'^ 

cm/s for T = 10 K and 7 = 5/3). Since shell 4 is 
10914 AU thick, sound wave will take 5.68 x 10^^ s 
{ti) to cross it. We assume that after every 
second, some percentage matter would move to 
shell 3. Similarly, we calculate the time scale for 
transferring some fraction of matter from shell 3 
to shell 2 (ta = 1.80 x 10 ^^ s) and shell 2 to shell 
1 {t 2 = 5.68 X 10^^ s). Due to outflow, some frac¬ 
tion of matter of shell 1 would contribute to shell 
4. Time scale for the outflow (ti) is chosen by con¬ 
sidering free fall time scale (\/37r/32Gp). Free fall 
time scale is inversely proportional to the square 
root of density. Here, we choose a typical dense 
cloud number density ( 10 "^ cm~^) for the calcula¬ 
tion of free fall time. So in our case we assume 
ti = 1.62 X 10^^ s. For simplicity, we assume that 
Pi = P3 = P2 = Pi = 20%. 


In Fig. 3, depth dependence of peak and final 
adenine abundances are shown for the non-mixing 
and mixing cases with a density profile of Fig. 2. 
Abundances are represented with respect to the 
total hydrogen nuclei in all forms. For each shell. 
We carry out our simulation up to 2 x 10® year. By 
the peak abundance, we mean the peak value ob¬ 
tained during the simulation regime and by the fi¬ 
nal abundance, we mean the value obtained at the 
end of simulation time scale. Since from Fig. 1 , 
it is clear that the pathways of Merz et ^ ( 2014h 
is dominating over all other existing pathways for 
the production of adenine in interstellar region, 
here, fo r Fig. 3, we cons idered only the path¬ 


ways of Merz et al. ( 2014ll . Clearly, due to the 


mixing of matter among various shells, final abun¬ 
dance of adenine is severely affected in the mixing 
model. Peak abundances of adenine is decreas¬ 
ing inside the cloud. Reason behind is that as we 
penetrate inside the cloud, density increases and 
this results in shorter depletion (to the next shell) 
time scale. Due to heavy depletion, all related 
species depleted from the gas phase. As a result, 
the peak value of adenine is reduced. However, 


the final value goes up and indeed, the final and 
the peak values become similar due to significant 
mixing. 


3.2. Results of Spectroscopic modeling 

3.2.1. Vibrational transitions 


HNCNH is one of the isomers of NH 2 CN which 
is a planar pentatonic molecule. It ha s nine fun¬ 
damental vibrations ( Rubalcava 1956ll . Similar 
to NH 2 CN, HNCNH also has nine fundamental 
vibration lines as shown in Table 3. In Fig. 
4a-c, infrared absorption spectra of HNCNH for 
the gas phase and grain phase are shown. We 
note that the gas phase HNCNH has its strongest 
feature at 2225.48 cm“^ and the second intense 
peak arises at 899.71 cm“^. In the grain phase 
2225.48 cm“^ peak is shifted to 2179.37 cm“^ 
and peak at 899.71 cm“^ is also shifted to 925.14 
cm“^ with much lower intensity. Beside these 
intense peaks, there are several peaks which are 
pronounced in our grain phase model (such as 
peaks at 902.61 cm“^ and 3564.13 cm“^). For 
the sake of better understanding, we placed all 
the infrared peak positions with their absorbance 
in Table 3. Moreover, a comparison with the 
other theoretical/experimental vibrational tran¬ 
sitions for HNCNH (Birk & Winnewisser, 1986; 
King & Strope, 1971) are also shown in Table 3. 
Some of our peak positions are very close to the 
results reported in the literature. 


In Table 3, we also note down the peak posi¬ 
tion of the infrared spectrum of C 3 NH in the gas 
as well as in the grain phase. We find that the 
most intense mode in the gas phase appears at 
2293.10 cm“^. This peak is shifted in the left side 
in the grain phase and appears at 2309.68 cm“^. 
Among the other peaks in the gas phase, 3718.05 
cm“^ and 463.33 cm“^ has the major contribu¬ 
tions. One strong peak at 3760.19 cm“^ appears 
in the grain phase. Fig. 4c shows infrared spectra 
of adenine in gas and grain phases. Gas phase in¬ 
frared spectrum of adenine contains strong peaks 
at 1634.40 cm“^, 1656.77 cm“^, 31.23 cm“^ and 
all these peaks are shifted to grain phase and 
appears at 1638.84 cm“^, 1621.56 cm“^, 141.19 
cm“^. Not only that, few new strong peaks ap¬ 
pears at 3594.23 cm“^, 3634 cm“^, 3718.16 cm“^ 
for grain phase adenine. These features are clear 
in the Fig. 4c. 
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3.2.2. Electronic transitions 


4. Concluding Remarks 


We continue our computation to obtain the 
spectral properties of C 3 NH and HNCNH in the 
electronic absorption mode. Electronic absorp¬ 
tion spectra of HNCNH, C 3 NH and adenine in 
gas as well in the grain phase are shown in Fig. 
5a-c. Corresponding electronic transitions, ab¬ 
sorbance and oscillator strengths are also summa¬ 
rized in Table 4. An electronic absorption spec¬ 
trum of HNCNH molecule in the gas phase is char¬ 
acterized by five intense peaks. These five transi¬ 
tions occurred at 219.55, 184.83, 156.97, 141.48, 
102 nm with major contributions from 1-Al —>■!- 
B2, 1-Al^l-Bl, 1-A1^1-B2, 1-A1^1-B2 and 1- 
Al—>■ 1-A2 electronic transitions. All the peaks 
are shifted in grain phase with the correspond¬ 
ing change in electronic transitions, absorbance 
and oscillator strength. For C 3 NH molecule, elec¬ 
tronic absorption spectra in gas as well in the grain 
phases are shown in Fig. 5b. This electronic spec¬ 
trum in gas phase contains peaks at 182.6, 126.7, 
104.04 nm due to 1-A' 1-A', 1-A' 1-A”, 

1 — A' —>■ 1 — A'. As in case of HNCNH, here 
also corresponding gas phase peaks shift in the 
grain phase. However, most interestingly, the gas 
phase as well as grain phase adenine contains only 
one strong peak in the electronic mode shown in 
the Fig. 5c. Gas phase strong peak appears 
at 256.12 nm with oscillator strength 0.0655 for 
1 — A —>■ 1 — A electronic transition. This peak is 
shifted to grain phase and appears at 252.11 nm 
with oscillator strength 0.1267 for the same elec¬ 
tronic transition. 

3.2.3. Rotational transitions 

In Table 5, we summarize our calculated rota¬ 
tional and distortional constants of the precursors 
HNCNH and C 3 NH of adenine and adenine itself, 
in gas phase. Calculated constants are corrected 
for each vibrational state as well as vibrationally 
averaged structures. Here we use B3LYP/cc- 
pVTZ level to perform these calculations in gas 
phases. In Table 5, calculated distortional con¬ 
stants correspond to I* representation with ‘A’ re¬ 
duction. Figure 6 shows the stick diagram of two 
new precursors of adenine. 


Interstellar chemistry to form complex bio¬ 
molecules began in this century and the first quan¬ 
titative computation of adenin e and other pre- 
biotic m olecules was presented by 
( 2000al fhll . A follow up study by 


Chakrabarti & Chakrabarti 


Maiumdar et al 


(j 2012 n showed that formation of adenine mainly 
follows radical-radical/radic al-molecular r eactio n 
pathways as proposed by ICupta et al.l ( 201 lh . 
More recently. iMerz et al.l ( 20141) performed retro 
synthetic analysis to find out a new mechanism 
for adenine formation in the gas-phase. They pro¬ 
posed that it could be formed from three already 
known interstellar molecules, namely, C 3 NH and 
the isomers, HNCNH and }l 2 NCN. Uniqueness of 
this pathway is that it does not involve HCN, wa¬ 
ter or ammonia during the production of adenine 
by any of the six intermediate steps. In this paper, 
we carried out a comparative study among various 
formation pathways o f adenine. We considered all 
the pat hw ays used in Chakrabarti fc Chakrabarti 
(l2000al lbl) : ICupta et ah ( 20 111) : Merz et al. ( 2014 ). 
Our chemical model suggests t hat the aden i ne for - 
mation pathways proposed bv I Merz et al. ( 2014 ) 
dominates over all other existing pathways. We 
studied the effects of mixing caused by both ad- 
vection and outflows. We find that for adenine, 
though the peak abundance is severely reduced, 
the final abundance becomes higher as compared 
to the non-mixing (static evolution) case. We 
explored the infrared, electronic, sub-millimeter 
spectroscopy of adenine along with its two new 
precursor molecules (HNCNH and C 3 NH). We 
show that adenine will have one single strong line 
either in the gas or in the grain phase vibrational 
spectra. A detailed chemical and spectroscopical 
information about adenine as revealed by our anal¬ 
ysis and its precursors would be extremely helpful 
for future survey of this species in the interstellar 
medium. In case the adenine is not directly de¬ 
tected, we believe that the detection of HNCNH, 
H 2 iVCA^ and C 3 NH would give adequate reason 
to believe that adenine must also be present. 
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Table 3: Vibrational frequencies of precursors HNCNH and C 3 NH of adenine in gas phase and H 2 O ice 


containing grains at B3LYP/6-311G++(d,p) level 


Species 

Peak positions 
(Gas phase) 
(W^avenumber in cm~^) 

Absorbance 

Peak positions 
(water ice) 

(Wavenumber in cm~^) 

Absorbance 

Experiment or other 
theoretical value 
(W^avenumber in cm~^) 


540.31 

79.926 

539.48 

141.348 

537^ 


540.65 

0.195 

539.60 

0.142 



718.17 

107.962 

719.18 

140.079 



897.11 

9.975 

902.61 

712.664 

890 ± lO'’, 886*^ 

HNCNH 

899.71 

438.296 

925.14 

18.842 



1287.51 

0.0593 

1284.96 

0.3371 

1285 ± 20 1275*’ 


2225.48 

713.448 

2179.37 

1290.062 

2104.7*’, 2097*“ 


3596.90 

143.912 

3564.13 

269.387 



3599.41 

26.106 

3569.93 

52.647 



180.31 

3.324 

154.87 

207.6288 



190.52 

1.491 

187.05 

38.291 



463.33 

409.108 

188.75 

2.815 



592.72 

2.186 

563.39 

8.437 


C 3 NH 

596.46 

62.223 

572.20 

2.130 

- 


965.33 

0.683 

953.06 

23.805 



1957.85 

27.241 

2028.00 

259.828 



2293.10 

1589.606 

2309.68 

2393.2106 



3718.05 

448.645 

3760.19 

1561.688 



^ Birk and Winnewisser, 1986 




^ King and Strope, 1971 




Table 4: Electronic transitions of precursors HNCNH and C 3 NH of adenine at B3LYP/6-311++g(d,p) level 


theory in gas phase and H 2 O ice containing grain phase 


Species 

Wavelength 

(gas phase) 
(in nm) 

Absorbance 

Oscillator strength 

Transitions 

W^ave length 

(H 2 O ice) 

(in nm) 

Absorbance 

Oscillator strength 

Transitions 


182.6 

47462.573 

0.9608 

l-A’^l-A’ 

191.57 

61654.567 

1.5193 

l-A-s-l-A 


126.7 

12259.477 

0.0073 

1-A’->1-A” 

143.78 

10621.323 

0.0276 

1-A-s-l-A 

C 3 NH 

104.04 

9869.950 

0.0216 

1-A’^-l-A’ 

119.07 

11042.926 

0.0593 

1-A-s-l-A 


- 

- 

- 

- 

104.3 

7174.6216 

0.0034 

1-A-s-l-A 


219.55 

672.854 

0.0166 

1-A1-S-1-B2 

215.94 

716.869 

0.0177 

1-A1-S-1-B2 


184.83 

1043.151 

0.0217 

1-Al-s-l-Bl 

159.3 

17358.816 

0.4126 

1-A1-S-1-B2 

HNCNH 

156.97 

10575.579 

0.2225 

1-A1-S-1-B2 

144.73 

22602.596 

0 

1-A1-S-1-A2 


141.48 

22671.816 

0.4362 

1-A1-S-1-B2 

123.29 

8869.585 

0.1003 

1-Al-s-l-Al 


102 

18576.3854 

0 

l-Al-i.l-A2 

100.94 

14075.884 

0.0338 

1-A1-S-1-B2 


- 

- 

- 

- 

87.24 

5237.201 

0.1028 

1-Al-s-l-Bl 


Table 5: Theoretical rotational parameters of adenine and its two precursors HNCNH and C 3 NH at 
B3LYP/cc-pVTZ level of theory 


Species 

Rotational 

Values 

Distortional 

Values 


constants 

in MHz 

constants 

in MHz 


A 

378941.67576 

Dj 

0.30806x10"*" 

HNCNH in gas phase 

B 

10550.84447 

Djk 

0.36739x10“ 


c 

10334.08855 

Dk 

0.11423x10“ 




di 

-0.89721x10"“ 




d2 

-0.12535x10"** 


A 

1431947.28270 

Dj 

0.48833x10"“ 

C 3 NH in gas phase 

B 

4694.90368 

Djk 

-0.19558x10“ 


c 

4680.33880 

Dk 

0.41486x10“ 




di 

-0.20704x10"“ 




d2 

0.29323x10"® 


A 

2387.05922 

Dj 

0.19731x10"** 

C 5 H 5 N 5 in gas phase 

B 

1575.31931 

Djk 

0.14017x10"“ 


c 

949.02096 

Dk 

0.10519x10"^ 




di 

-0.10350x10"** 




d2 

-0.49530x10"“ 
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Fig. 3. -Radial distribution of adenine when mixing and no 

mixing models are considered. Due to heavy mixing the final 
and peak abundances (dashed lines) become similar. 




Fig. 4. - Infrared spectra of HNCNH, C 3 NH and adenine 

in gas and grain (water ice) phases. 
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Fig. 6 . - Rotational stick diagrams of HNCNH, C 3 NH and 

adenine. 


Fig. 5. - Electronic absorption spectra of HNCNH, C 3 NH 

and adenine in gas and grain (water ice) phases. 
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